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ABSTRACT 


This thesis presents an interactive naval radar model 
which computes radar detection in the presence of land 


masses, using a parametric terrain descripticn. 
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I. INTRODUCTION 
This «hesis presénts a novel naval radar model which 
computes radar detection in the presence of land masses. The 

model is an interactive computer program which accepts 


scenarios and radar parameters from the user and displays a 


map of the area indicating where targets can and cannot be 


(D 
ry 


@evectved. Th esulting map can be displayed at the user's 
computer terminal or printed offline. 
Major capabilities of the model are: 
f. pavasgetric terrain description 
2. Wser friendly input and output 
3. beach return masking 
WY. radar shadowing 
5. side lobe masking 
[r= program 2S Written in FORTRAN IV H (extended), to be 
axecuted on an IBM 3033 using an IBM 3278.2 video computer 
terexneal. & data file which contains a parametric terrain 
description nust be prepared pefore using this program, Ddut 


all other required input is prompted at the terminal. 





ieee AOTMATION 

The Warfare Environment Simuiator (WES) used by Command, 
Communications, and Ccntrol personnel is a large scale con- 
puterized naval wargame. The land dispiayed by the systen 
Seeseamoc adxtrcct the Eadar detection probabilities. This is 
Ome Of the artificielities in the game. 

During Exercise BRIGHT HCRIZON '81 the tactic of con- 
cealing small vessels in fjords and among islands to prevent 
redar detection until a surprise attack could be made was 
Gemonstrated to be effective. This technique, while it 
Sommompe anticipated, can not currentiy be modeled in WES. 

Wezth Atlantic Treaty Organization naval warfare Sgenac- 
ios Virtually always include the proximity of land. United 
States Navy scenarios are primarily open ocean. It seems 
necessery to incorporate NATO problems into USN procedures. 
With this in mind, a method of encoding terrain data was 


combined with a radar model to produce a method of producing 


25.2 Onew2 en tue Concealed areas 
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III. MODEL 


by 


a. WARFARE ENVIRONMENT SIMULATOR 


W=S is a computer program that Simulates a naval 
warfare environment in a large computer system, enabling 
personnel to engag2 in realistic wartime scenarios without 
she expense of actually losing ships and lives { Ref. 1]. 

2. Radar Equation 

The WES model uses the standard radar equation modi- 

fied =o handle ducting and clutter effects. Equation (1) is 


used to calculate signal excess. 


SE=P++2G+2W+TCS-4DR-B-NF-L-C (1) 


62) 
UJ 
iH 


Signal excess, dB 


yg 
qt 
HW 


peax transmitted power, dB//w 
G = antenna gain, dB 
“q = wavelength, dB//cn 


TCS= target cross section, dB//n 


Pe e=ecuctihg cactcr 

R= target tange, d3B//na 

B = seceiver noise figure, dB 
He = noise factor, dB 
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Gy = Gaanter factor, dsb 


a 
il 


system loss factor, dB 

system loss factor includes antenna pattern loss and 
atmospheric absorption loss. Antenna pattern ioss takes 
into account «he change in illuminating energy lavels 
prough= apou*« by the lobe shape of the radar main beam pat- 
tern. The internal system losses are assigned to be 1.5 dB, 
a typical value for mcst radars. Atmospheric absorption loss 
2s added to this. Atmcspheric absorption loss is frequency 
Gevendent, and is assigned as 1 dB around 300 4Hz and 3.5 dB 
az 5000 Mdz, act target ranges typical to naval radars. 

The clutter Lactor describes the losses due to sea 
Cl@rcer. Tt is taken to be 10log(Hw), where Hw is the signi- 
Ficaae wave height in feet. 

Preweveceeng) Cackeos is used to describe the effects 
Oe Glek=-a@ propadation, especiaiiy surface evaporation ducts. 
An approximate fit of the IREPS model of a leaky wave guide 
Ss eed. * A strength of zero indicates no ducting, anda 


Bace2 hor. 260 amitaticn is evoxed. If a duct is present, 


wo 


fhis restriction is lifted. <A strength of five corresponds 


Mmenmeegneted Refractive Effects Prediction oY Seen 
eloped by Naval Ocean Systems Center and inplemented on 
fewest ~Packard Model 9845 desktop calculator. It 4s a 
board Siiviron#ental data processor and display systen 
gned +o aid in the assessment of the impact of lower 
Spheric refraction effects on Naval EM Systems. 


MM fulM ct Qa 
t(D tr oD 
BW td a 
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to a "perfect", lossless duct with free space propagation 


Within its poundaries. Values for D are selected from Table 


Table ls" Ducting Factcr 


cee OC titi‘ i = 2 
oie Ss 0 ee | a ha 


The noise factor is used to include jamming. It is 
met wWSeed in this ecdel. 

See subsection 3 below for the method of processing 
Signal excess. The method used in determining the radar tar- 
Met CEOSS section of ships is discussed in Chapter III, sec- 
mon Cc, swbsection 3. - 

For «his model an additional term is included. The 
Pobstacie gain" is calculated and subtracted from the right 
Rand side of the equation.= 

3. Signal Excess model 

Signal excess is converted to probability of detec- 

tion on the basis of false alarm number and pulse integra- 


Zion as di ad below. 


in 


Cc 


i) 
ib 


= 


bay 


j~ 


alse alarm number is the numoer of radar signal 
Butses that are received before a non-Signal soise pulse is 


2For obstacle gain see chapter III, section c, subsec- 








received. WES, as well as this mcedel, offers the use of four 


Gitferent values. 


—_ 


The number of pulses illuminating the target is 
evaluated by use of equation (2). 


60 (PRR) (Bw) 


Np= ----------- 
(S®) (ARR) 0 


wheres; 
ND = number of pulses 
PRR = oulse repetition rate, pulses per second 
BY = horizontal beam width, in degrees 
S# = angular width of swepth sector, in degrees 
ARR = scan rate, in scans per minute 


The distribution or signal excess can be closely 
approximated by a normal distribution with a mean defined by 
the number of pulses integrated, probability cf false alara, 
and a standard deviaticn of 7dB. A large takle contains the 
mean ofr the signal excess distribution (au) in terms of nun- 
Der Of pulses integrated and probability of false alarn. 


Conversion cf computed signal excess (SE) to prob- 


3 


ee cy Oc detection is accomplished in WES as follows. 


(a 


~ 


The signal excess is normalized by equation (3). 


= (se-mu) /sd (3) 


N 


The z is used to enter a table of cumulative standardized 
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normal values to get the probability of detection which 
tseused gm the WES simulation. 

This model runs the process in reverse: given the 
desired probability of detection, the number of pulses, and 
an assumed standard deviation the standardized normal prob- 
ability distribution is used to find a corresponding sigzal 
excess, which can «hen be compared to the realized signal 
excess tO determine if the target will be detected within 
*he constraints given. 

4. Conditions and Assugptions 
A number of constants are established in the WES 


model, as well as in this modified version. Many of then 


ela] discussed in the sections in which they are singu- 


larly avpropriate, but some of the more general ones will be 
Mentioned here. 


A 4/3 earth radius is used in calculating the dis- 
Bane=) sO the radar horizon to compensate for standard atmnos- 
pitee#c refraction. 

Receiver noise figure is taken as a constant value 
of 5.5 dB, a value which is primarily thermal noise under 
Seesearda conditions. 

The return frem land is not directiy modeied. The 


geometrical cross section cf a portion of parametrize 


As. 





3 


terzain is assumed to be the radar cross section. Data is 
Unavailable on the reflectivity of varying types of surface 
20 compute the true radar cross section given the gecmetric 


cross section, so a reflectivity of 1 is assumed. A slope 


O 
tty 
fv 
ct 
ss 


east .02 1s reguired for beach return, when checking 
from the seaward side, and a slore of at least .3 is 
reguired for return ir the sidelobes. In ¢€ither case, the 
land must be at least as high as the target vessel. These 
values were determined by experimenting with the geometrical 
cross section of a terrain sample. The steeper slope 
tequizred for the sidelobe masking is caused by the lower 


gain present in the sidelobe. bs 


Be SMO LER-@ILLS MODEL 


eu 
Q 


7. Basic System 
In 1979 Josef Smoier wrote his Master's thesis on an 
“Operational Lanchester-Type Model of Small Unit Land Con- 
bee [ Ret. 2). It is a time sequenced, deterministic, batta- 
lion-level, force-cn-force model implemented on a digital 


NS 


Pa aon 1d 


tu 


cf nod 


a 


AG. ve 


(f 
(D 


cOmpucer. I+ contained a aetho 
oped by Christopher James Needles in March, 1976 for his 
Master's Thesis "Parameterization of Terrain in Army Combat 
Analysis" [Ref. 3]. In September 1980 Glenn MW. Mills 


attempted to overcome shortfalls in the original model, and 





added several enrichments to provide added user flexibility 


tty 


(Ref. 4]. A user's manual is provided for this model on a 
pPezt@anent disk in the W. R&R. Church Computer Center. These 
Fortran programs conduct their Army battles over land, using 
Needisas' terrain model and Professor Hartman's elevation and 
line-of-sight routines.3 These routines were used to model 
the land in this naval model. Only minor modifications were 
required, partially tc adjust to the change in scale and 
Pact@eailly to remove scme strictly land effects. 


2s 


}+ 


et ba. 


1) 


ode 
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in vot6 
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hristoph¢er James Needles presented and 
evaluated a methodology of parameterizing terrain for use in 
Wad cé@bat amalysis {Ref. 3}. This was a shitt from the 
theme -pornal method of Gigitizirg data compiled from terrain 
and interpolating. He described a metnod by which terrain 
could be created mathematically by using a modified bivari- 
ate normal probability density function. 

The common fern of the bivariate normal density 


function is given by equation (4). 


=—_— <i ae ee eee ee ee ee ee ee ee ee ee ee ee ee ee 


We 





: | 
| 0 0 (9) | uct nu'0" 
(1) | f---=1)) “eparagende-Lae-) SxXepde | gy (eu lee ayy 


CO 
cd 





The equation has been modified by making the normal- 
izing constant equal the the naximum elevation of the 
desized terrain. The resulting equation has sufficient par- 
ameterzs to model many different types of "hills", and by 
carefully combining these hills different types of terrain 
Gan be aodeled. 

This method of terrain modeling is currently being 
Usecery she STAR model weitten in SIMSCRIPT I1.5, but the 
Fortran version of the method has also been used, in the 
Smoler-Mills nodel vréeviously described. ¢ 

Soemmegqure (1) £or an illustration on hoW the hills 


are construczed. 


The puppose of the function ELEY is to compute ter- 
Eeaimeel=Yaleion for given xX, Y coordinates. 

The elvation routine was initially provided to Josef 
Smoler by Professor James K. Hartman for use in his land 


combat mod2l, and subsequently by Mills in his modification 


thy 


aa | 
ry 
iD 
th 


cx Smolerts nodel ico |ts 


( 
rt, 
wr 


The point defined Dy the Coordinates is exanined to 


determine if any of the "hills" are present, and if so which 


see@e acicn of Taczical Miternative Response, a 
ground-air combat model developed at the Naval Postgraduate 
SGmool. 
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| 
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: tae,Shape of =he hilis from above is : 
| modified by the correlation coefficient | 
| 
| Le TOP VIEW ) 
me yl! ! 
| ep | 
, : 
| | X | 
: | 
| (i, | 
) SIDE VIEW 
| 3 | 
| | | 
| 
| , | 
| | 
| ) 
| | | 
: | 
l 
| END «VESA | 
: | 
| 
| ( 
| | 
| | 
| 
| 
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mmm mmm a ae ee eee eee eee ee ceil eee 


Egqunmees:. Terrain Structure 


ones. Only those which are effective are evaluated, mininiz- 


ing <=n= computer time required fof each point. Given the 
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Bngex of the “effective hill", the appropriate paramenters 
are selected and applied +o the previously defined nodified 
ivariate normal distribution to derive the terrain eleva- 

meome at thatepoint. 

12a gags] cond elevacion routine, ELEVG, the gradient 
components of the equation are computed as well as the ele- 
vation, allowing the slope to be determined easily. 

The procedure to add tree height to the top of the 
hills to increase elevation has been removed for this radar 

4. Line of Sight Routine 

The line-of-sight routine was also developed by 
Professor Hartman [Ref. 5]. Its development and use follow 
closely those of the elevation routines. 

Line of sight (LOS) is a purely geometric computa- 
tion assuming verfect visibility. The result of the LOS conm- 
Dutation was initially the percent of the vertical height of 
the target visable to the searcher, and the percent visible 
Of Wee Ssarcher by the target. In addition, the cange and 
Gummer On of che point waich causes the target to be totatly 
Gbeteueted 1S now returned. 

The basic procedure is: "Find the lowest sight line 


#EOG@ 4He searcher over the terrain. Extend this line to the 


ZA 





target's positicn and compare its extrapolated height to the 
~arget's elevation. Thus ccmpute the percent visible." 

In aS much as the searcher and the target are both 
Ships, the searcher is assumed to be at sea lavel. If duct- 
ing is present, the target is also assumed to be at sea 
tewel. If ducting is not present, an additional routine con- 
Fuces the amounts by which both the obstructing terrain and 
the target are lowered by the earth's curvature. 


As in ELEV, ferest calculations have been removed. 


. AUGMENTATION 
1. Znatroduction 

The STAR LOS and ELEV reutines have keen combined 
Witn =he WES radar equation to get a routine to calculate 
Geaee=2On behind and next to land. 

A separate routine was written to calculate beach 
SeteceS, Side lebe effects, and a special routine to fro- 
vid2 output on an IBM 3278 computer video terminal. 

After accepting preliminary terrain, tadar, target, 
and anvironnental data th2 program samplés sites a desig- 
nated distance apart from a specified reference point to 
create a gwatrix containing either elevation or detection 
d@ta, to be cxamined at the terminal, with an option to 


Dapee it On a line printer. 
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riguze 2: Flow Chart 


Figure (2) presents a flow chart of the overall 


prograa. 


ate 


tH 


n utines 


to 


we 
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(iD 
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Ollowing subroutines were developed to input 


tj} 


+he ianformation required by the programs. 
ap 46 fSEbaln Data 


tems are required to 


ts 


ferent 


tt 


A number cog di 
Gefine the terrain. The elevaticn of the reference plane 
@ise be entered, along with the total number of “hills”. For 
each hill, the X and Y coordinates, its maximum elevation, 


*he standard deviation in both the K and Y direction, and 
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Mee ™eea@rclatlion between them fMust be entered for each hill. 


This defines the hili, in as much as these are the apprepri- 


ct 


ate parameters in the bivariate normal Gistributioa 
described in the secticn on the terrain model.5 

The data for the terrain is assumed to be con- 
tained in a data file accessable py the virtual machine. Two 
Powe2nes adapted from ones writted py Mark L. Yount perforn 


the necessary file definitions, after requesting filename, 


Filetype, and filemode (filemode defaults to a) at the 


Information for up to 100 hills covering over an 

area 200nm by 200nm may be entered this way. 
b. Radar Data 

Radar information may be entered in either of 
BmOrmclis: First, it may be stored in a data file, and loaded 
pimemtcn tO the Mephod utilized to load terrain data. Sec~ 
ond, it may be entered at the terminal, in response to ques- 
tions on the screen. Table (2) lists the information 


bequ- fed. 


r 3YN equation used in the terrain model, see chap- 
peseceicn 8, Subsection 2. 
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fee 2s ekeguired Radar Bata 


NAME OF RADAR SYSTEM 

PEAK TRANSMITTED POWER (gaTTs) 

RECEIVER IF 3ANEDWIDTH (#42) 

PULSE REPITION BATE (PPS) 

HORIZONTAL BEAM WIDTH (DEGREES) 

VERTICAL 3EAM WIDTH (DEGREES) 

ANGULAR WIDTH OF SWEPTH SECTOR (DEGREES) 
SCAN RAT2 (3?) 

PULSE LENGTH (4ICROSECONDS) 

FREQUENCY (#2272) 

ANTENNA HEIGHT ( 227) 

S-E= COORDINATES OF RADAR (NM) 

NUMBER OP LOBES (UP TC 5) 

STRONGEST # LOBES (DEGREES CW FROM MAIN LOBE) 


(note: the "#" above represents the number of lobes desired) 


Co Gasgqet Data 
The target data is entered in the same subrout- 


ime. Table (3) lists the required information. 


Table 3: Target Data 


TDENEIFIER FOR TARGET 
VemeICAL SIZ OF TARGET Ae Pe) 
TeBO@2T MAXENUM DISPLACEMENT (ALTLOTONS) 


Additional information is required for the pro- 


gram. It consists of environmental and probabilistic data. 


Taoae 4: Bevrronneneal and Probablliistic Data 


DUCTING STRENGTH (INTEGER-O FOR NO DUCT TO 5S FOR PERFECT DUCT) 
PROBABILITY OF FALSE ALARM 

RFQUIRED PROBABILITY OF DETECTION FOR SINGLE SWEEP 
SIGNIFICANT WAVE HEIGHT (FEET) 


Table (8) lists these addit Honecl requirements. 
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A covy of anitial conditions is then available 
to be printed on the line printer upon reguest. Figure (3) 


is) am example of this listing. 
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Figure 3 





Finally, information relevant to the actual 
plotting is requested. Simply enter the number of nautical 
miles between sample points desired and a reference point on 
che overall map and the routine commences its calculations. 

3. Computing Target Cross Section 
The Target Characterization Sranch of the Naval 
Research Laboratory has nade carefully controlled measure- 
ments or the radar cross sections of a number of naval 


Ships. A simple empirical expression has been obtained (see 


equation (5)) that expresses the cross section as a function 


orf the displacement and the frequency. 
+o — tat 
e@s=52f°d (5) 
Where: 
tcs = Target cross section is square meters 


$ radar frequency in megahertz 


d displacement in kilotons (full load) 
This expressicn provides reasonable approximations 


for the microwave band common to most naval radars, and over 


he displacement range of 2000 to 17000 tons [Ref. 6]. 
4, Beach Return Sesking 


= 





Each segment of land (elevation greater than 0.5 
meters) is checked to determine if it is a beach. This is. 


done in subroutine BEACH by checking every adjacent sample 
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pose@= af the grid and determining if any of them is ocean 
feome-MaesoOn less than 0.5 meters). If this is the case, it is 
assumed that a potential target can get next to it. Next 


Subroutine BRIN estimates the geometric cross section of a 


0) 
(D 
Q 
ct 
1? 
O 
J 
oO 
rh 


beach. This is done by first determining the area 
within one resolution cell. The size of the cell is deter- 
mined by taking the area between two circles at the given 
range separated by the range resclution and between two 
Tadial lines separated by the angular resolution. Subroutine 
ELEVG is then called ~o determine the slcpe of the section 
of beach at the point of ccncern as seen by the radar. The 
effective area of the resolution cell is determined by tak- 
amg the sine of this Slope and multiplying it by the resolu- 
taem cell total area. See Figure (4) ror an illustration. 

Ii the effective cross section determined in the 


reater than that of the target vessel, it 


iD 


above @anner is 


‘Q 


is assumed that ary return from the beach would mask the 


return from the target. 


LAND is 2 Subroutine which checks the main lobe and 


given major side lobes for the presence of land. A proce- 


dure siniliar +o BEACH is used +o determine the proximity of 
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Seen io besear- wcneexea DY USIng ELEVG =o check for 
jand az the same range as the point to be evaluated in the 


direction of a major lobe when the target is centered in the 
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WR any of the inspected locations indicates that masking 
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fakes place. 
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ling the line-of-sight routine, the progran 
ee, COntac® sn sone of three ways: 


a. Case 1 


mi ees SeCOmOpGeructson. The obstruction gain is 


mee, Case 2 
were Ts Onlymcare al ObStruction. The targe= is 
snierfiace between shadow and illumination. The obsta- 
Mecs see at the log of an eighth, in accordance with 


Bemep=-Cal Daserference technique. [Ref. 7] 


Ge Cas2 3 
Pao aadeo 2st locallyropszructed, in which case 
routine OBGAIN is called to determine obstacle gain. 


OBGAIN calculates the "obstacle gain" created by 
lo waves being ditfracted over the hill tops. A 
edg2" obstacle is assumed, and the calculations are 


ed using Picquenard's radio wave propagation equa- 


Mee. 9). %ITe the hill ts not the highest of the 


the nill, and the target it is assumed that the tar- 


again at +h2 interface and tae same value mentioned 
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+urned, ctherwise the calculations continue. 
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Ih any case, the obstacle gain is treated as an 
attenuation and subtracted from the other values in the 
equation used to calculate signal excess. 
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PHA is a Wersion of THAP modified to allow the use 


rt, 


or spécial symbols. TMAP has been submitted for the NONIMSL 
LIBRARY at NPS. The routine interactively takes an array and 
Pempamees a Contour gap Zor the vicinity of a given point. 
Tepmomeeoie 2S ON a teeminal or on the orinter, at the option 
oz the user, interactively taking input from the keyboard. 
Passed arquments: A is the array to be plotted 

Nand M ar? its dimension 

SFNM, is the numbez Of nautical miles 

between points. 

Other arguments are requested at the terminal. the 


Beogmam 2s a nodifZed version of Professor Gilles Cantin‘s. 


Srogram CTRMAP version of OCT 1, 1969, adapted for terminal 
Figure (5) is given on the screen immediately before 


commencing the plots, either initially at the terminal or 


printed at the top of the paper printout. 
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Symbology for Plot 


5 


Figure 


Note: On Figure(6), the "3" representing the radar's position is off the Screen. 








Figure (6) is an example of what appears on the ter- 
minal. The paper listing gives the entire map area, with 


symbology as before. 


| 
| PLOT CENTERED ON 45.00NM SOUTH AND 45.00NM ZAST | 
| sie teak ke A 
| whe BBR A 
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‘#a8CD F G De 
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Figure 6; Screen Display 


In the screen display given, the points are one nautical 
mile apart, with the radar located approximately 20nm to the 
mBOrthwest of the center of the island. The letters roughly 
centered in the screen correspond to terrain elevation, as 
SGedecamin Begquz> (2). The '*®* on the left edge correspond to 
a steep beach, where che beach return is stronger than the 
reflecticn from the target. The '®* to the right define the 
area shadowed py the island. The clear area around it is sea 


level. 
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This reutine all 


owsS the model to be rerun after 


changing a select number cf the parameters. Figure {7) is 


the screen display which is the menu of allowable changes. 


Simply enter the number 
nd then the new value. 
mew varameters and then 


elewakeion and de=t=etion 


corresponding to the desired change 
Tae ZOuUtine will offer te prine the 
Gouse > uct another Nadtrexecontaining 


data for examination. 
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IV. ZUTURE 


=. PesolBke UTILIZATION CE MODEL 

This model has applications in planning and training for 
Naval Operations. 

mee training purposes, this model could be incorporated 
into WES or other sea combat models to better reflect detec- 
tion in areas where land is present. For such uses, the 
input and display routines could be triamed, and only the 


point in question would be analyzed, not an entire map 


For planning purposes, ship courses could be selected to 
provide maximum radar protection from land shadowing, or the 
inverse, to allow minimum terrain interference for station- 
t@g Fadar guard shios, could be found after parameterizing 


the applicable areas. 


de POPURE SNRANCEMENTS 

Darabolic cylinders could be utilized to model the 
=efFraction effects instead of knite-edges to improve the 
nodel somewhat (Ref. 9]. 

Pa=  inpm@e and ou€put routines could be modified to take 
better advantage of facilities available (plotting 
terminals, for instance). 


a7 








Bpemetieod co find tne beach, instead of identifying if a 
Pel kewrear point is within a given distance of the sea, 
would enarle better ccastai coverage. 

Radar return is often observed behind obstacles under 
conditions in which simple knife~edge refraction cannot be 
responsible. The causes for these effects shouid be found 
ag@ incorperated into the model. 


For use in other models, different subroutines may be 


th 


"4 


adapted for the particular nodel being constructed. Some 
Subroutines will have to be entirely rewritten, especially 
those involving beach return which use mapped data already 


calculated, to handle single-pcint problems. 
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mee foloving is arn ertine listing of the radar nodel. 
Be bas been specially written to be run on the 18M 3033 on a 
IgM 3278.2 video terminal, the facilities available at the 
Naval Postgraduate School. Appendix B contains and explana- 


ron Of Gach of the subroutines. 
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The follewing subroutines are use NDAs apeejgect: 


&. SBACH 

Thas Subroutine dstenmines if a point of land is adjacent to 
*he sea. 

Pe 2nd g@6 avrewere K and YY coordinates in array A 

A is the map array, dimensioned N by M& 


KOAST is the returned value 


B. BRIN 
This routine calculates the geemetric cross section of an 


StecmoOewatiaetne size Of a resolution cell which may contain 


Kigt1 ate the positicn of the radar 
K2yt2 age the positicn or the beach 
RRES is the range resolution 

Me25 is the angular resolution 

TMAC is the terrain elevation 


ECS is the effective radar cross section of the beach 
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This subroutine computes the elevation reduction caused by 
the earth's curvature present at 4 given range. 

BER is the effective earth radius in meters 

RANGE is the distance traveled along the earth 

DROP is the loss in elevation resulting from the earth's 


curvature 


This routine determines the elevation at a given location. 
apr eare the coordindtes of the point 


TMAC is the returned elevation 


tJ 


PWhoisempcutine is Simaliar co ELEV, but it also calculates the 
gradient components. 


GX AND GY are the calculated components of the gradient. 


E. “GeTss 

This routine processes radar data to compute the signal 
PaGees cremmemed to ke radar from 4 contact under the speci- 
ia conditions. 


SE is the returned signal excess 
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detection data, 


sie. 


INPOT 


= 
- 


out 


(t 


am> reads in the redar, target, environmental, 
either from the terminal or from a data 
he variables are listed in order of appearance 


ah) 


Sher subroutine) 


x 


PTR 


py 


DRS 


ARES 


Sa 


ARR 


NP 


TG 


sa 


name of the radar systen 
peak transmitted power 
receiver bandwidth 

Dulse repetition rate 
horizontal beam width 
vertical beam width 
angular resolution 
antenna gain 
angular width of swept sector 
scan rate 

pulse length 

Frequency 

radar antenna altitude 


- ~~ 
— af io 


~adar antenna al ude 


Geordinmates cE radar 
wavelenath 


nuaber of pulses illuminating a target per sweep 


identifier for target 
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ECS 


i Sam 


ANF 


REIN 


RM AX 


RHORIZ 


tt 


Mertical size of target 

vertical size of target 

*arget maximum displacement 

target radar cross section 

ducting strength code letter 

Dpeovadwtsty O©zutaciure code letter 

reguired probability of detection for a single 
sweep 

Signal excess required to give a desired PDET 
Significant wave heigat 

Cite eEs Factor 

system loss factor 

receiver noise figure 

Minimum range 

range resolution 

NMaxinum unambigious range 


Bacar HOri.zon 
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Er. @8OVER 


This routine is only used internally by LOS 


J.  DASED 
This subroutine determines if a point would be masked by 
return from land in either the main beam adjacent to the 


Given point, or by land at the appropriate point ina side 


MASK is the returned value, one if there is masking, zero 


otherwise. 


KK. LOS 

[ees TCOuUtIne calculates the line-of-sight in terms of a 

Buete=loOn Bistple for observer-target pairs. (the variables 

are listed in the ordér of appearance in the subroutine) 

mm Ye (82,Y8) are the X,¥ coordinates on the field for Aa 
and 3. 

PRACA (TMACB) are terrain elevation for A (8B) 

PHeCkh (7TMICB) clevation of A (B). 

SIZBA (SIZEB) vertical dimension of A (B). 

BATOS (LEEOA) ardicator variable €écr LOS calis. 

MepPevsmcne herman. Of che obstzucting hiil (if any). 

RANGE is the range from the radar to the opstructing hill 


(1f-a ny). 
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VISPRA (VISFR3) fraction of SIZEA(SIZEB) which can be seen 


Divi BA). 


L. O38GAIN 
This routine calculates the “obstacle gain created by the 


tadio waves being defracted over hill tops. 


H1 elevation of the radar 

HO elevation of the hill 

5 2levation of the target 

D1 range from the radar to the hill 
D2 range from the hill to the target 


RLAMDA wavelength 


OBGN returned obstacle gain 


%. SPARS 
Tres SOUcsee aS ™by Mark Yount. It splits a double-precision 


string down into its separate words. 


Ba RANGE 
Mires LroOutame finds che range between (X1,¥1) and (%¥2,Y¥2) in 
@a plane. 


RANGE 1s returned. 


78 





Ge “Ra AP 

This routine prepares contour maps with special symbols to 
be displayed either on the terminal, in response to loca- 
ions keyed in, or on a line printer. RMAP is composed of 
three subroutines: PREMAP, PSTMAP, and PRTIMAP. 

A is the array to be plotted, dimensioned N by M 
SPNM is the number of nautical miles between sample 


Per»mts on the grid. 


Pen Soieoce 
This routine defines the data files so that the program may 


access «hen. Weemten Wath assistance from Mark L. Yount. 


oO. SrPvp 
This routine reads in the terrain cata from a predesignated 
Gata file. (the variables are listed in the order they 
appear in the subroutine) 

MiELiS teal nu@ber cf hills 

BASE Surface reference plane 

XC(I) and YC(I) are the x,y coordinates of the center of 
tiem as 1 

PRAK(I) is the maximum elevation of hill I. 

SX(I) and SY({I) are the standard deviations in X and Y of 


Mzl1 fi. 
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RHO (TI) 
Boe 


NHL 


NHTOT 


BPSii 


GOr=sclation beeween K andeyY for hill [I 

DOinter On where .o start treading on LISTH 

azray holding how many hills aze in corresponding 
Grid Square 

total number of hills in LISTH 


a list of which hills affect which gridsquare 
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